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Abstract 



K^ , We have investigated the effect of the isospin i, J^ = 0^ resonant state iCo(1950) 

L •• I on the decays D^ -^ K^rj and D^ — > K^rj' as a function of the branching ratio sum 

cn ■ r = fir(K^(1950) -^ K^r]) + 5r(K^(1950) -^ K^r]') and coupling constants gx^ROr,^ 

9K*K°ri'- We have used a factorized input for D^ — > i('g(1950) weak transition through 
Q>^ , a Tri^ loop. We estimated both on- and off-shell contributions from the loop. Our 

0^ \ calculation shows that the off-shell effects are significant. For r > 30% a fit to the decay 

,^ ' amplitude A[D^ — > K^rj') was possible, but the amplitude A{D^ — > K^rj) remained at 

♦^1 its factorized value. For small values of r, r < 18%, we were able to fit A{D^ — > K^i]), 

Qh! and despite the fact that A{D^ — > K^rj') could be raised by almost 100% over its 

factorized value, it still falls short of its experimental value. A simultaneous fit to 
both amplitudes A{D^ — > K^rj') and A{D^ -^ K^rj) was not possible. We have also 
determined the strong phase of the resonant amplitudes for both decays. 

rS 
^ . (PACS numbers: 13.25.Ft, 13.25.-k, 14.40.Lb) 



I. INTRODUCTION 

A common method of evaluating the matrix elements of two-body hadronic decays of 
heavy mesons, B and D, is based on the factorization approximation [|I| which utilizes 
model form factors. However, this approximation has had only limited success in describing 
two-body hadronic decays of the D meson |0, |^, ^, |§, §]. In particular, the factorization 
approximation not only underestimates the decay rates for D° — >■ K^t] and D'^ — ^ K^rj', it 
generates r{D^ -^ K^rj) > r{D^ -^ K^rj') in contradiction with the experiment 0. In an 
attempt to remedy this discrepancy, Ref. P studied the above decays in the factorization 
approximation but included the annihilation term. They found that unlikely large form 
factors for K — > ^7(77') transitions were required in order to bridge the gap between theory 
and experiment. Ref. on the other hand introduced nonfactorized contributions and used 
a flavor-S'[/(3) parametrization for the nonfactorized matrix elements to fit the data. Their 
conclusions imply a large value for the hair-pin amplitude. 

Hadronic decays of mesons are complicated by the presence of final-state strong interac- 
tions {FSI) between hadrons in the final state. The importance of FSI in hadronic decays 
of D meson has been known for a long time; its role was emphasized by several authors [§| 
almost twenty years ago. The long-range FSI generates phases in the decay amplitudes 
and the most dramatic effect of FSI is induced by the interference between different isospin 
amplitudes which depends on the phase difference [|I], |], |I^ between different isospin ampli- 



tudes. In decays with a single isospin final states, as in D^ — * K^rjlr]'), isospin interference 
effects are absent. However, FSI also leads to a change in the magnitude of the decay am- 
plitude, and not simply a rotation in the complex amplitude plane. Hence we expect FSI 
to effect the decay rates in single-isospin channel decays too. 

In this paper we have considered the FSI effect of 7^0(1950) resonance on D^ — > K^t] 
and D^ -^ K^rj' decays. The mechanism we are proposing for resonant FSI is as follows: 
D^ — >■ K^ri{ri') are color-suppressed decays in the factorization approximation ( see Fig. |l|). 
The resonance ii'Q(1950) has a substantial branching ratio (~ 50%) to Ktt mode, leaving 
room for its coupling to Krj{ri') channels. We propose that the effect of 7^0(1950) on D^ -^ 
K^ri{ri') could be estimated via the Feynman diagram shown in Fig. |^ where in the loop we 
include both K~t:^ (color-favored decay) and K^ti^ (color-suppressed decay) states. Such 
mechanism has been invoked in Ref. |]ll] in D^ —* K^K^ decay. However, in contrast to 
Ref. |TT]], where only the on-shell loop contribution is retained. We evaluate both on- and 
off-shell loop contributions in Fig. ^. 

Recently, Gronau ll^ has also discussed the role of resonant FSI on D decays. We relegate 



a discussion of these works [|TT|, [T2[ , and their relationship to ours, to the last section of this 
paper. 

In Sec. II we detail the model and the method of calculation. The results are presented 
in Sec. HI, and Sec. IV deals with the discussion. 

II. METHOD OF CALCULATION 

A. Calculation without final-state interactions 

The decays D^ -^ K^rj, K^rj' are Cabibbo-favored and are induced by the effective weak 
Hamiltonian which can be reduced to the following color-favored {CF) and color-suppressed 



(CS) forms § (G = ^K.K*d): 



HcF = G[ai{ud){sc) + C20s], 
Hcs = G[a2{uc){sd) + ciOs], 



fll 



where Vqq> are the CKM matrix elements. The brackets (ud) etc. represent {V - A) color- 
singlet Dirac bilinears. Og and Og are products of color-octet currents: Os = ^ J2l=i {uX"'d){sX"'c) 
and Og = I J2a=i (mA^c) {sX°'d) where A*^ are the Gell-mann matrices. The parameters ai and 
02, related to the Wilson coefficients, are chosen to be ai = 1.26 ±0.04 and 02 = —0.51 ±0.05 
0. Note that these values of Oi and 02 are achieved in A^ ^ 00 limit, where N is the number 
of colors. 

In the factorization approximation one neglects the contribution from Os and Og, and 
the matrix element of the first term is written as a product of two current matrix elements. 
Since we are effectively working with N y^ 3, one could argue that the nonfactorization 
arising from Og and Og is being included since working with A^ 7^ 3 is equivalent to working 
with A^ = 3 but including nonfactorization effects. The quark diagram which contributes to 
D^ — i> K^ri{ri') decays is shown in Fig. |l|. The physical particles 77 and rj' are mixtures of the 
flavor-singlet tjq = -^^luu + dd + ss > and the flavor-octet rjs = -j^\uu + dd — 2ss > with a 
mixing angle 6p = —20° ITEI: 



rj = rjs cos 9p — 770 sin 9p 
r\ = r]8 sin 6p + ijq cos 6p. 



(2) 
(3) 



The factorized amplitude for the decay D^ 
'factorized') 



K^f]{f]') is given by (superscript / refers to 



Af{D^ -^ K^T]) = Ga2 {K° \ sd \ 0) [t] \ uc \ D°; 
In calculating the decay amplitude in Eq. (^), we use the following deflnitions, 

{K{p)\isd),\0) = -tfKp, 



{'n{p,^)\{uc)^\D{pD)) 



PD+Pr, 



mjj-m; 



F^\q') 



(4) 
(5) 
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k,F,^\q 
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where q = pd — p^^ \s the momentum transfer, fx is the decay constant of the K meson, 
Pi ^{q^)-i (^ = O5 1) ^^6 invariant form factors deflned in [|14|. The factorized amplitudes for 
D° -> K^t]{t]') are 



Af{D^ 
Af{D' 









sine'iml- ml) F^\ 



cos^'(m|) 



riJ- U \''^K. 



K) } 



(7) 



where 9' is given by 



sin^' 



v^ 



cos 9 1 



sin 9 J- 



COS 6*' = \ -cosOp + \ -sinOp. 



(8) 



The corresponding decay rates are given by 



r^(D° ^ K'^viv')) =1 A^ 



ircm 



2 ' 
D 



(9) 



where | p | is the center of mass momentum and mp, is the mass of the decaying D meson. 



B. Calculation with resonant final-state interactions 



The resonances contributing to D^ -^ {K^Vi K^v') niust have isospin and spin-parity 
assignment I{J^) = \{0'^)- There are two particles, ^0(1950) and ^9(1430), with such 
properties [0: mi = tuk* = 1945 ± 10 ± 20MeV, Ti = 201 ± 34 ± 79MeV and ms = tuk* = 
1429 ± 4 ± 5MeV, Ts = 287 ± 10 ± 21MeV. Ahhough ms + ^ is much smaller than m^ 
and ^^0(1430) decays almost exclusively to iiK channel, we cannot prejudice its effect to be 
insignificant in D^ —>■ {K^ri,K^r]') [15|. However, in this paper we consider the contribution 
from ir*( 1950) only. 

The contribution of ^q(1950) to D^ -^ K^ri{r]') is represented by Feynman diagram 
shown in Fig. |^, where in the loop we consider both the color-favored state i^~7r+ and the 
color-suppressed state K^tt'^. Evaluation of the diagram in Fig. ^ gives the following ampli- 
tudes ( superscript r and subscript — h, 00 refer to 'resonant' and K^ir^, K^tt^ intermediate 
state, respectively) 



a: 



;{l/2) 



d^k- 



V_i'>{k' 



{2tiYJ {k^ - ml){{w - kf 

I-+ X A!_*r^ 
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A strong 



(10) 
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^str^ng ^ ^ 



nip — rap) + iVmp 



/, 



(11) 
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A strong 



(12) 
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4 strong 
^00 



5'007 — 2 2 \ , -n ^' (1^) 

[nip — nijf) + li mp 

where g^j^ and Qqq are the couplings of K~7t'^ and A'^vr^ states to -^"0(1950) and / is the 
coupling of K^rjijj') state to ^q (19^0) 5 ^^R is the resonance mass, F its width, w is the four 
momentum of the decaying particle [w = {mp, 0) in the CM.) and k is the loop momentum 
to be integrated over; V_+ {k'^) and Vqo (^^) are the vertex functions in isospin 1/2 state. 



They are related to the amphtudes A{D^ -^ K^n^) and A{D^ -^ K^tt^) and are evaluated 
in the following. 

Although A!l+ and Aqq get contributions from color-favored and color-suppressed inter- 
mediate states respectively, they are not independent but are related by isospin and SU{3) 
symmetry as the following analysis elucidates. As the resonance i^o(1950) has isospin 1/2 
and as strong interactions conserve isospin, only the isospin component A^^^"^^ in the fol- 
lowing will contribute to Vi^ (^^) and Vijf (^^)- The isospin decomposition of the decay 
amphtudes. 



AiD'^^K-n^) = J1a(3/2) + j|a(i/2) 



A{D^^K%^) = W^A(=^/2)- J1a(i/2), (14) 

allows us to calculate A^^^"^^ which is needed as the input. For the color-favored intermediate 
state K~7i~^ we get from Eq. (11^ 



v'S\e) = J^A(V2) 



Ia{D' -. K-n^) - ^^-A{D' -. K%') } , (15) 

while for the color-suppressed intermediate state, we get 

I^AiD"^ -^ K-1T+) - ^-A{D^ -^ K\'') } . (16) 

Therefore, isospin symmetry offers the following relation between the two vertex functions 

Vif)(fc2)andVir(^') 

V^X"\k') = ~V2Vg/'\e). (17) 

Let us now turn to the coupling of Kq{1950) to mesons. Using a S'f/(3)-invariant strong 
Hamiltonian which couples a scalar octet S to two pseudoscalars octets P, 



HsPP = '-Tr{{P^,P^}S), (18) 



we obtain 



3 9 /.ON 

9 = 9K*K-tt+ = 9-+, 9K*K0n0 =900 = -7^, 9K*K0rjs = /^- {^^) 

From equations ([Tl|), (|13]) and (|19D the SU{3) symmetry gives the following relation between 

,1 J tj- J istronq i tstronq 

the strong amphtudes A_^ '^ and Aqq ": 



Hence, using equations (p!OD, (p!2D, (p!7|) and (^) the resonant amplitudes A'^__^_ and Aqq are 
related by 

A% = 2AI, (21) 

and the total resonant amplitude A^ is 

A = A__^_ + Aqq 

= l^L^- (22) 

To proceed further with the calculation models and approximations are used. 

For an explicit calculation of the amplitude ^1+ in Eq. ([I0|) we need to calculate the 

integral /_+ for which we require the momentum dependence of the vertex function V_+ (k"^)- 
For this we assume the form dictated by the factorization assumption but with both tt and 
K not necessarily on their mass-shells (see Ref. [|16]), 



Af {D' ^ K\p)^\k)) = -rG^fK{ml-e)F^\p% (23) 



From Eqs. (^) and (|23|) we then obtain 



;(i/2)/,2n _ „-rr,/^ I ,/2^ . . 2 /„ u\2M:^DK. 




V.Tik') = -tG\l-{\l-aMmj,-iw-ky)Frik') 

~^fKiml-e)Ftiiw-kr)y (24) 

where p and k are four-momenta of K and vr mesons in the loop, respectively. They are 
related by momentum conservation at the vertex, w = p + k. Since in our calculation, the 
intermediate particles {K, it) are allowed to be off-shell, we have to assume a behavior of the 
form factors in the vertex function, Eq. (plD, as the particles go off-shell. Form factors with 



a dipole dependence: F{k\) = I ^2^^ ) (where A is fixed by experiment) have been used 



^x 



in the past [T^ to describe off-shellness of intermediate state particles. In this work we have 



used the following phenomenological form factor 



^0 (^ ) - ,. ^2..^ (^_fc)2-™,2. . ■ (25) 



(1-|^)(1 



A2A-^ A2 



-x 



In principle, one could use two different mass scales, Ai and A2. We chose to work with a 
single mass scale for simplicity^. Beside having a dipole dependence, this form factor satisfies 
the following limit: As K goes on-shell ( p^ = (ly — /c)^ — > {m^ = '^'k) the form factor in 
Eq. (pSf ) reduces to the usual monopole form, 

i^o^^(fc^) = ^P^, (26) 



^The mass scale A in the factor Fi^^ is different from the A in the form factor Ff^'^ . 



with the pole mass A given in Ref. []I4| . 

The integral / ^. is obviously complex, the real part arising from the region where both 

TT and K are off-shell and the imaginary part coming from the region where both tt and K 
are on-shell. The details of the calculation are provided in the Appendix. The values of 
the mass scale A are given in the next section. The final result is (to make the factorized 
amplitude real we drop a common factor of i from the amplitudes in Eqs. (^ and (p^)) 

/_, = 10-G^Fo--(0)|8.735 + 1.769^;^ + .(8.328 + 2.211^;^ 

= \I.+ \exp{t6i) GeV (27) 



Now, model predictions |]18| as well as experiments |]19| give a ratio pdkJ) ~ 1- Consequently 



the phase 6i ~ 45° and it is insensitive to form factor-models. The magnitude of / ^. depends 

on F(f^(0). We use the value F^^{0) = 0.76 |TB|, |g to obtain 

/_+ = 0.723 X 10^^exp(i45°) GeV. (28) 

The resonant amplitude A^ also depends on the strengths and the signs of strong coupling 
constants / and g which we determine as follows. The decay rate of a scalar particle decaying 
into two pseudoscalars is given by 

r(^ ^ pp) = ^14^. (29) 

STrrrig 

For ^q(1950), we have the following measured branching ratio [^ 

s = Br{K* -^ Kit) = 52 ± 14%. (30) 

Using Eqs. (0), (^) and the central value of s in Eq. (|30|), we get g = ±2.707 GeV. Since 
the determination of / = gK*K°r](r]') = 9ri{ri') is complicated hj rj — r]' mixing and the fact that 
no measurements are available for the branching ratios Br{KQ -^ Kr]) and Br{KQ -^ Ki]'), 
we provide some details of how we calculate (7^ and gr^'. We include r] — rj' mixing in the 
strong decay of the resonance as follows, 

gr^ = gs cos 9p - go sin 9p 

gri' = gs sin Op + go cos 9p, (31) 

where the octet coupling g^ is determined from Eqs. ( [T^ ) and (|TP|) to be gg = gK*K'^-qg, = 
—g/V^ and the unknown singlet coupling is go = gK*K^r)o- We treat the following unmeasured 
branching ratio sum as a variable, 

r = 5r(ir*(1950) -> Kr]) + 5r(ii'*(1950) -> Kt]') 

where p and p' are center of mass momenta of the final state particle in the decays Kq -^ Krj 
and Kq —>■ Kf]', respectively. Since gs is knownQ from flavor-S'f/(3) symmetry, we use Eqs. 

^The results of this paper were obtained using the positive value g — 2.707 GeV. We learn nothing new 
if the negative value of g is used. See Discussion. 
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(pT|) and (^) to solve for the singlet coupling go{r) in terms of r. We obtain two solutions 



for go, denoted by gQ{r),i = 1,2; we then substitute 5'o(^) in Eq. (^) and get two sets 
of solutions {gl^{r),gl^,{r)),i = 1,2. Their dependence on r is shown in Figs. |^ and ^. In 
order for the strong coupling constants grj{r) , grj' [r) to be real we find that we must have 
r > 5%. We also, have the constraint r + s < 100%, which restricts the allowed range for r 
to: (5% < r < 52%). Using Eqs. (|10D and (|19D with g = 2.707GeV, the amplitude /li*;°"9 is 
found to be 

A'T' = -5.45^;(^,)(r)exp(^52°). (33) 

Finally the total resonant amplitude A"^ is (in the following calculation all the amplitudes 
are scaled by a factor of 10~^), 

= -3.9%;(^,)(r)exp(297°) GeV, (34) 

and the total resonant amplitude A^ , 

A' = -A' 
2 -+ 

= -5.91^;(^,)(r)exp(z97°) GeV. (35) 

Note that the amplitude A^^^""'^ defined in Eq. (|10D is complex; however, the phase of A^ is 
not the phase of A^^^"""^. The total amplitude for D^ -^ K^r]{ri') is the coherent sum 

A = A^ + A'. (36) 

III. RESULTS 

For numerical calculations we used the following parameters, 

Vcs = 0.974, Vud = 0.975, U = 0.131 GeV, Jk = 0.160 GeV, (37) 

and the pole mass A, 

m(sc(0+)) = 2.6 GeV, m(rfc(0+)) = 2.47 GeV (38) 

The values of the form factors are presented in Table I. 

A. Factorized Amplitude 

For the purpose of comparison with experiment we have calculated the factorized am- 
plitude A'l'lD^ — > Kr]{ri')) using the following models for the form factors: i) Bauer, Stech 
and Wirbel (BSW) model |]14|, where an infinite momentum frame is used to calculate the 
form factors at q^ = 0, and a monopole form for g^ dependence is assumed to extrapolate 
all the form factors to the desired value of g^; ii) Casalbuoni, Deandrea, Di Bartolomeo, 
Feruglio, Gatto and Nardulli (CDDGFN) model |^, where the form factors are evalu- 



ated at g^ = in an effective Lagrangian satisfying heavy quark spin-flavor symmetry in 
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which hght vector particles are introduced as gauge particles in a broken chiral symme- 
try. A monopole form is used for the q^ dependence. The experimental inputs for this 
model are from the semileptonic decay D — > vr/z/, and we have used the value of the form 



factors F]"^(0) = F^^{Q) = 0.83 ± 0.08 |21[ extracted from data, and decay constants 



/d, = 213ln ± 11, /d = 194l}o ± 10 MeVlgg in calculating the weak coupling constants of 
the model at g^ = pOl , which are subsequently used in evaluating the required form factors 



; vi) Isgur, Scora, Grinstein and Wise (ISGW) model [^, where a non-relativistic quark 
model is used to calculate the form factors at zero recoil and an exponential g^ dependence, 
based on a potential-model calculation of the meson wave function, is use to extrapolate 
them to the desired g^; iv) Lubicz, Martinelli, McCarthy and Sachrajda (LMMS) model |P4|| , 



where the form factors are obtained from lattice calculation of D meson semileptonic decays. 

The factorized amplitudes with model form factors and the experimentally determined 
amplitudes are presented in Table II. The prediction for D^ — > Kr] amplitude is too low in 
every case; an enhancement of a factor of 1.5 to 3, depending on the model, is needed to 
match the experimental amplitude. For D^ -^ Krj' , the situation is worse. 



B. Resonant Amplitude 

In the following we list the amplitude represented by the diagram of Fig. ^ separately 
for the cases where the loop particles are on-shell and off-shell. 



1. On-shell contribution: 

The contribution to the resonant amplitude A^ from on-shell loop particles is obtained 
by taking the imaginary part of the integral J_+. We get from eq. (p2D, ( pH]) and p3|), 

A'ion-sheU) = hlm{I^+) x Al*^^ 

= -4.18^;(^,)(r)ea;p(zl42°) GeV. (39) 

2. Off-shell contribution: 

For the loop particles off-shell the resonant amplitude A"^ is obtained by taking the real 
part of the integral /_+. We get, 

A"- {off -shell) = ^Real{I.+) x A!*™"^ 

= -4.184(^,)(r)ea;p(i52°) GeV. (40) 

Note that the on-shell contribution, Eq. (|39|) , is of the same size as the off-shell, Eq. (|40|) , 
but advanced by a phase of 90° which comes from the factor i. The on-shell and off-shell, 
amplitudes have the same magnitude but different phases, therefore including off-shell effect 



will modify both the amplitude and the phase of the resonant amplitude. The total resonant 
amplitude A^ is given by, 

A' = A' {on-shell) + A" {off -shell) . (41) 

Finally, the total decay amplitude is obtained by adding the factorized amplitude A^ to the 
resonant amplitude A^, 

A = Af- 5.91^;(^,) {r)exp{i97°) GeV. (42) 

Plots of the magnitude |y4| = |y4-^ + A*"! as a function of r are given in Figs. |^ and ^ 
for the decays D -^ Krj and D —>■ Krj', respectively. In these figures we have used 
Af{D° -^ Krj) = 8.37 GeV and A^(D° -^ Krj) = 4.5 GeV as predicted in BSW model 
only. We discuss the results in the next section. 



3. Strong phases 

The phase of the total resonant amplitude and the on-shell amplitude can be determined 
from Eqs. ( P3| ) and (|5UD, respectively. The sign of the coupling constant, g^^ and g^^, is 
important; a change in sign introduces a phase difference of 180° which leads to a different 
pattern of interference (constructive or destructive) between the factorized amplitude A^ 
and the resonant amplitude A^. We use the graphs of Figs. ^ and ^ to read the signs of gri{r) 
and grj' (r) and then determine the strong phase of the resonant amplitudes. The phases are 
summarized in Table III. 



IV. DISCUSSION 

The factorization prediction for the amplitude A{D^ —>■ K^if) is too low compared to the 
experimental data and the situation is even worse for A{D^ -^ K^r]') (see Table II). 

Amplitudes with the iclusion of resonant FSI (Eq. ^) are plotted in Figs. ^ and ^ for 
A{D^ — > K^if) and A{D^ —>■ K'^r]'), respectively. From Fig. ^ we notice that for solution 
i = 1, A{D^ -^ K^T]) ~ Af{D^ -^ K^Tj) over most of the range of r, except where r is low, 
r < 14%, when A^lD'^ -^ K^rf) gets a small enhancement over its factorized value, but stays 
below the experimental value. On the other hand. Fig. |] shows that for solution i = 1, 
A{D^ —>■ K^rj') rises with r and fits the experiment in the range 30% < r < 42%. Although, 
for lower values of r A{D^ -^ K^r]') is underestimated, it still gets large enhancement over 
its factorized value. Solution i = 1 does not accomodate a simultaneous fit to A{D^ -^ K^rf) 
and A{D^ -^ K°t]'). 

As for solution z = 2, we notice from Figs. |^ and ^ that both A{D^ —>■ K^rj) and 
A{D^ — > K^rj') rise with r. In particular, a fit to A{D^ -^ K^t]) is secured for 6% < r < 20%. 
Despite the fact that in this range of r A{D^ -^ K^r]') could be raised by almost 100% over 
its factorized value, it still remains well below its measured value. Again a simultaneous fit 
to A{D° -^ K°T]) and A{D° -^ K^t]') is eluded. 
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Recently Dai et al. |]ri| have used the same mechanism, but kept only the on-shell loop 
contribution to estimate the effect of resonant FSI in the decay D^ — > K^K^ . Hence, con- 
trary to our case, the strong phase is solely determined by the resonant propagator. 



1^ uses current algebra 
K*(1430)7r 



with 
measured 



In a recent paper Gronau |T2[ has calculated the contribution of i^o(1430) in the direct 
channel (annihilation topology, as in our case) to D^ —>■ K^n^ decay in a model- independent 
way and found it to be a substantial fraction (~ 20%) of the total amplitude (which is 
largely isospin 1/2). He argues that the effect of i^Q(1950) on D^ -^ K~tt^ decay could even 
be larger. The difference between his approach and ours (apart from the fact that we are 
dealing with different D^ decay modes) is that we have used the factorized input for the weak 
transition D^ -^ ^^0(1950) through a t^K loop, while Ref. 
smoothness assumption, to relate D^ -^ ^q(1430) vertex to D 
rate. Thus while the D^ — >■ ^q(1950) vertex of Ref. [1^ is assumed to be real, ours is 
complex. If our resonant contribution is equated to the W-exchange amplitude of Ref. [|I2|, 
then clearly the resonant contribution has a phase around 90° (modulo tt) relative to the 
tree amplitude (see Eq. (^)). 

In summary we find that the resonant FSI due to i^Q(1950) in the direct channel effects 
A{D -^ K^rj) and A{D —>■ K^rj') significantly. However, a simultaneous fit to both decay 
amplitudes is not possible. 

A final comment: We tried the same calculation with a negative sign for g, g = —2.707 GeV. 
Apart from leading to a phase shift of 180° in the strong phase 6r, it did not change our con- 
clusions. It was still impossible to fit both A{D — > K^t]) and A{D — * K^t]') simultaneously. 

Appendix 

Because of Eq. (|2^) , the calculation of the resonant amplitude A"^ reduces to that of 

A'^_^ = I |_y4!_!|°"^. Using the vertex function given in Eq. ( p^ in Eqs. (|3^ ) and (|I2|) 

we get the following integral 



-iG 2 r d'^k 

^2^V si {k^ - ml){{w - kf ~ m\) 



X 



^\a,Uml -{w- kf)F,^^{k') - ^Mml - k')F,^-{{w - kf) } , (43) 



where d'^k = dkod^h, k"^ = k^ — |kp, and w = {mu,0). Using form factors as explained in 
the text, the integral / |_ reduces to the following form. 




V6 
where the integrals Jj are given by 



^ { an |AV.Fo^^(0) Uml - m^)/f ^ - Jf ^ 



-a,l^A%^^0)\{ml-ml)I^--I,^- 



(44) 



h 



(A2-m2)(A2-m2^) 



(-^^AA — IaK — In A + Ittk) 



(45) 
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D° 




K' 



7], 1] 



Figure 1: Quark diagram contributing to the factorized amplitude A^ for D^ 



0^ 1^0^/ 



X^ri, K\ 




V, V 



Figure 2: Resonant contribution to D^ -^ K^rj, K^rj'. The vertical dashed line represents 
the cut when the particles in the loop are on-shell. The thick line represents the resonance 
-^"0(1950), and the shaded circle represents the weak vertex V{k'^). 



h 



-1 



-1 



-I TV A + 



1 



(A^ — m 



Iak + 



(A^ — ml 
1 



--^^AA, 



27-* AX 



(A2 - ml 



fT^AA- 



The integrals li = Jj with A having the appropriate mass given in Eq. 

IxY has the generic form 



iXY 



dknd^k 



(A;2 — mx){{w — ky — my)' 



(46) 
(47) 

38]) . The integral 
(48) 



First we integrate over dk^ in the complex plane with a contour closed in the lower half plane. 
We have used a mass scale Eq. (|38D A > m£,o + m,K, therefore the integration over d^h is 



well defined for all the integrals except I^^k which has a pole contributing to the imaginary 
part of /_+. The above calculation leads to the following result, 



'- {8.735F(f^(0) + 1.769F(f"(0) 
-«(8.328Fo^^(0) + 2.211Fo^"(0))} x lO^^G GeV. 



(49) 
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Figure 3: Plot for the two sets of solutions, (?* 



1,2, as a function of the branching 



ratio sum r = Br{K^{1950) -^ K^rj) + 5r(ir;;(1950) -^ K^i). The thick and light parts of 
the curve correspond to gHr) and giij-), respectively. 




Figure 4: Plot for the two sets of solutions, gL{r),i = 1,2, as a function of the branching 
ratio sum r = Br{K^{l950) -^ K^rj) + Br{Kl{l9bQi) -> K^r]'). The thick and light parts of 
the curve correspond to gl^iij-) and g^i{r), respectively. 



Table 1: Model predictions for the form factors 
for the processes D^ — > K^t], K^t]' , Kir. 



(0) 






^0 {^k)i 







(0) and the ratio 





BSWI 


CDDFGN 


ISGW 


LMMS 


Fo'^'iml) 


0.710 


0.313 


0.638 


0.344 


Fo^'^'iO 


0.683 


- 


0.937 


0.240 


Fo^^(O) 


0.762 


0.699 


0.769 


0.63 


Fo^'^(O) 


0.692 


0.83 


0.510 


0.58 


Fjf-(O) 


0.91 


1.19 


0.66 


0.92 
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factorization 



Figure 5: Plot of the magnitude of the total ainphtude \A\ = \A-^ + A^\ in Eq. ( ^21) for D^ — >■ 
K^f] as a function of the branching ratio sum r = Br{KQ{1950) -^ K^t]) ± i?r(i^Q(1950) -^ 
K^Tj'). The thick and hght parts of the curve correspond to the solution g]^ and g"^, respec- 
tively. The shaded region represents the experimental value of the amplitude, the horizontal 
line represents the factorized amplitude A^ . The values on the y-axis must be multiplied by 
a factor of 10~^ GeV to get the absolute magnitude of the decay amplitude. 



25 



20 



10 



Experiment 




factorization 



Figure 6: Plot of the magnitude of the total amplitude |y4| = |y4-^±74^| in Eq. (^2p for D^ — >■ 
K^r]' as a function of the branching ratio sum r = i?r (7^0(1950) — > K^if) ± i?r(i^Q(1950) — »• 
K^r}'). The thick and light parts of the curve correspond to the solution g]^, and g"^,, respec- 
tively. The shaded region represents the experimental value of the amplitude, the horizontal 
line represents the factorized amplitude A^ . The values on the y-axis must be multiplied by 
a factor of 10~^ GeV to get the absolute magnitude of the decay amplitude. 



Table 2: Model predictions for the factorized amplitude A^ for the process D^ — > K^r]{r]') 
These values must be multiplied by a factor of 10~^G'e\^. 





BSWI 


CDDFGN 


ISGW 


LMM 


Expt. [0] 


At\D^ - 
Af{D^ - 




8.37 
4.50 


3.70 


7.53 
6.18 


4.06 
1.58 


11.3 ±0.8 
20.54 ±1.54 
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Table 3: Phases of the total resonant and on-shell amplitudes defined in Eqs. (^) and (|39|), 
respectively, as functions of r for the processes D^ — > K^ri{r]'). 6r is the phase of A^ (Eq. 
(PSP) and Son-shell is the phase of A^' (on-shell) (Eq. (|5D|)). 



Solution 


Decay 


r 6r Son-shell 


i = I 


^0 — , ;^0^/ 


(5 < r < 35)% 97° 142° 
(35<r<48)% (97 ±180)° (142 ± 180)° 
(5<r<48)% (97 ±180)° (142 ± 180)° 


i = 2 




(5 < r < 48)% 97° 142° 
(5 < r < 48)% 97° 142° 
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